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FOREWORD 

u*JupP 

The  present  study  is  part  of  a  program  of  "Theoretical  Research  on 
Combustion  and  Wave  Propagation  in  Metal  Dust-Oxidant  Mixtures "ybeing  con 
ducted  by  the  Division  of  Engineering,  Brown  University,  under  United 
States  Army  Ordnance  Contract  Wo.  DA-19-020-0RD-4761.  The  work  was  admin 
istered  by  the  Pyrotechnics  Laboratory  of  the  Picattiny  Arsenal,  Dover, 


N.  J 


SUMMARY 

A  macroscopic  theory  of  two-phase  flow  with  mass,  momentum  and  energy 
exchange  is  discussed.  The  theory  is  applied  to  the  study  of  systems  which 
depart  only  slightly  from  local  thermodynamic  equilibrium.  An  example 
of  wave  propagation  in  a  two-phase  medium  with  viscous  and  thermal 
relaxation  is  calculated.  It  confirms ' some  of  the  predictions  of  ref.  2, 
where  a  thermodynamic  analysis  was  carried  out. 


1.  Introduction  The  analysis  presented  here  was  undertaken  in  an  attempt 
to  understand  better  the  nature  of  propagation  of  waves  and  combustion  fronts 
in  a  pyrotechnic  mixture.  The  mixture  consists  usually  of  two  components, 
both  in  granular  form:  a  light  metal,  such  as  aluminum  or  magnesium,  and 
an  oxidizer,  such  as  potassium  perchlorate,.  The  granules,  have  an  average 
diameter  of  15  microns  which  is  usually  much  smaller  than  the  smallest 
macroscopic  scale  of  interest,  (see  ref.  l)  For  example,  in  a  particular 
problem,  the  smallest  macroscopic  scale  of  interest,  denoted  by  j[  ,  may 
be  a  fraction  of  a  millimeter  or  larger.  Furthermore,  there  are  usually 
so  many  granules  in  a  volume  Jl  that  the  mixture  may  be  considered 
macros copically  as  a  continuum.  The  implication  of  a  continuum  description 
of  the  granular  mixture  will  be  discussed  more  fully  in  the  next  section. 

For  the  moment  we  shall  accept  this  view. 

The  mechanical  and  thermodynamic  properties  of  such  a  continuum  de¬ 
pend  on  the  microscopic  structure  of  the  mixture;  in  particular,  on  the 
microscopic  state  of  the  aggregate  of  the  granules.  There  are  two  extreme 
cases  which  may  be  discerned:  In  the  first  case,  the  mean  distance  between 
neighboring  granules  is  large  compared  with  the  mean  diameter  of  the 
granules;  in  the  second,  the  two  lengths  are  of  the  same  order  of  magnitude. 
In  the  former  case,  the  great  majority  of  the  granules  are  not  in  contact 
with  another  granule.  These  granules  are  merely  particles  suspended  in  a 
perfect  gas.  Such  a  mixture  behaves  essentially  like  a  "Dusty  Gas".  In 
the  second  case,  most  of  the  granules  are  in  direct  contact  with  their 
neighboring  granules.  As  an  extreme  example,  the  granules  are  so  tightly 
packed  together  that  they  lock  into  one  another.  In  this  extreme  case, 
the  mixture  behaves  essentially  like  a  solid,  or  better,  a  "porous  solid" 
with  the  intergranular  space  filled  with  a  perfect  gas.  In  the  inter¬ 
mediary  range  between  the  "dusty  gas"  and  the  "porous  solid",  the  mixture 
assumes  a  mixed  behavior.  Thus,  for  example,  the  mixture  exhibits 
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properties  which  are  characteristic  of  a  solid  undergoing  plastic  defor¬ 
mation.  The  thermodynamics  of  a  dusty  gas  and  a  porous  solid  are  given  in 
refs.  2  and  3j  that  of  a  loosely  packed  granular  mixture  is  given  in  ref.  4. 

In  a  thermodynamic  analysis,  it  is  assumed  that  the  relative  motion 
between  the  granules  and  the  surrounding  gas  is  negligible  only.  Further¬ 
more,  we  assume  that  the  granules  and  the  surrounding  gas  are  either  in 
thermal  equilibrium  with  one  another  at  all  times  or  are  thermally  insulated 
from  each  other.  In  spite  of  these  rather  severe  restrictions,  the  result 
of  a  purely  thermodynamic  analysis  often  leads  to  important  conclusions. 

An  example  of  such  an  application  is  given  in  ref.  2. 

In  what  follows  we  shall  be  interested  in  problems  in  which  the 
granules  and  the  gas  between  them  are  not  in  thermodynamic  equilibrium. 
Between  these  two  components,  there  will  in  general  be  an  exchange  in  mass, 
momentum  and  energy.  We  propose  to  examine  these  non-equilibrium  phenomena 
within  the  framework  of  a  continuum  theory. 

Earlier  contributions  to  this  subject  were  mainly  due  to  F.  A.  Williams 
who  developed  a  statistical  theory  of  spray  and  applied  it  to  many  combustion 
problems,  (ref.  5*  6,  7)  We  adopt  here  an  alternate  approach  —  the 
macroscopic  approach.  We  shall  make  use  of  the  method  of  irreversible 
thermodynamics  to  establish  certain  phenomenological  laws  describing  the 
interaction  between  the  flow  and  the  particles  suspended  in  it.  The 
limitations  of  these  laws  are  discussed.  An  application  of  the  macroscopic 
theory  to  wave  phenomena  in  such  a  medium  is  included. 

2.  Some  Macroscopic  Properties  of  the  Mixture  We  shall  now  examine  the 
conditions  under  which  a  granular  mixture  may  be  considered  as  a  continuum 
on  the  macroscopic  scale  and  we  shall  discuss  some  of  the  properties  of 
such  a  continuum.  This  discussion  can  perhaps  be  best  illustrated  by  con¬ 
sidering  a  special  instance. 

Suppose  that  the  average  diameter  of  the  granules  is  of  the  order  of 


10  microns,  i.e.,  10”^  cm,  and  the  average  distance  between  the  neighboring 
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granules  is  of  the  order  of  100  microns,  or  10  cm.  It  follows  that  in 
1  cubic  millimeter  of  the  mixture  there  are  about  1,000  granules.  If  the 
smallest  macroscopic  length  of  interest  is  greater  than  a  half  millimeter, 
to  be  abbreviated  by  l(:0,5wm)  ,  then  all  linear  dimensions  equal  to  or 

less  than  X  and  all  areas  or  volumes  equal  to  or  less  than  and 
respectively  may  be  taken  as  macro scopically  zero.  Under  such  conditions, 
it  may  be  possible  to  define  the  various  macroscopic  density  variables  at 
a  point  in  the  mixture  considered  as  a  continuum.  Thus,  for  example,  the 
mass  density  of  the  granules  at  a  point  may  be  defined  as  follows.  Let  SV 
be  a  macroscopic  volume  containing  the  point  and  S M  be  the  total  mass  of 
the  granules  in  fV  .  The  ratio  SM/SV  as  S  £  (which  is  regarded  as 
macros copically  zero)  represents  the  mass  density^ of  the  granules  at  a 
point  of  the  mixture  considered  as  a  continuum.  Likewise,  the  stress 
vector  at  a  point  in  the  mixture  may  be  referred  to  the  force  acting  on  an 
element  <TA  as  etc.  On  the  other  hand,  to  define  a  mean  tem¬ 

perature  and  velocity  for  the  granules  and  the  gas  between  the  granules, 
we  must  exclude  very  rapid  changes.  An  estimate  of  the  smallest  macroscopic 
time  scale  allowable  for  defining  such  quantities  can  be  made  as  follows. 
First,  let  us  consider  a  single  granule  of  10  micron  diameter.  Assuming  a 

q 

thermal  diffusivity  of  0.2  cm  y  sec  for  the  granule,  it  is  seen  that  the 

1(P}  /0.2  seconds 

or  5  microseconds.  Hence  if  the  smallest  macroscopic  time  of  interest  is 
greater  than  5  microseconds,  we  may  assign  a  single  temperature  to  each 
granule.  Since  all  the  granules  in  a  volume  i>3  undergo  almost  identical 
changes,  we  may  in  fact  ascribe  a  mean  temperature  for  all  the  granules  in 
the  volume.  Similarly  a  mean  velocity  of  the  granules  may  be  defined  at 
each  point  of  the  mixture  considered  as  a  continuum.  Next,  let  us  examine 
the  fluid  surrounding  the  granules.  Since  the  average  distance  between 
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neighboring  granules  is  of  the  order  of  100  microns  or  10  cm.  and  the 
thermal  diffusivity  of  gas  is  of  the  same  order  of  magnitude  as  that  in 
aluminum,  namely,  0.2  cm  /sec,  therefore  the  thermal  relaxation  time  in 
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the  gas  between  the  granules  is  of  the  order  of  • 5x10  sec.  or  half  a 
millisecond.  Hence,  we  may  define  a  mean  temperature  for  the  gas  surround¬ 
ing  the  granules  when  the  smallest  macroscopic  time  of  interest  is  at  least 
1  millisecond.  Since  the  Prandtl  number  of  the  gas  is  of  the  order  of 
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unity,  the  viscous  relaxation  time  is  also  about  5x10  seconds  and  we 
may  ascribe  a  single  velocity  for  the  fluid  inside  a  volume  . 

It  follows  from  the  above  considerations  that  if  the  smallest  mac¬ 
roscopic  length  of  interest  is  large  compared  with  the  average  distance 
between  the  granules  and  if  the  smallest  macroscopic  time  of  interest  is 
long  compared  with  the  thermal  relaxation  time  of  the  fluid  between  the 
granules  then  a  macroscopic  continuum  description  of  the  granular  mixture 
is  feasible. 

There  is  one  further  point  to  be  clarified.  In  the  special  example 
discussed,  nothing  is  said  about  the  time  required  for  the  granules  and  the 
gas  surrounding  them  to  attain  the  same  temperature  and  velocity.  We  have 
only  discussed  the  time  required  for  defining  a  mean  temperature  and  velocity 
for  each  of  the  two  components  of  the  mixture.  An  estimation  of  the  time 
required  for  the  granules  and  the  surrounding  gas  to  attain  the  same  tem¬ 
perature  can  be  readily  made.  First,  we  observe  that  the  rate  of  heat 
transfer  from  the  gas  to  a  granule  is  of  the  order  of 


k.-AI_,.TT  D1 
3“D/l 


(2.1) 


where  K^is  the  conductivity  in  the  gas,  AT  is  the  temperature  difference 
between  the  granules  and  the  surrounding  gas  and  D  is  the  diameter  of 

the  granule.  .AT  is  thus  a  measure  of  the  average  temperature  gradient 

D/2 

at  the  surface  of  the  granule.  On  the  other  hand,  the  total  change  in 
heat  content  in  the  granule  as  a  result  of  heat  transfer  is 
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(2.2) 


j'rlTj]r-cPAT 

where  the  term  in  the  bracket  is  clearly  the  volume  of  the  granule  while 

and  Cy  are  'respectively  the  density  and  specific  heat  of  the  granule. 

The  time  required  for  the  establishment  of  thermal  equilibrium  between  the 

two  components  is  obtained  by  dividing  (2.2)  by  (2.1),  namely,  II  D1 

3  k  ^ai 

If  c*  denotes  the  thermal  diffusivity  in  the  granule  (^  d,  -  0-2  cw\a/sec^ 

This  timers  JL_._h _ t  P1  where  k  is  the  conductivity  of  the  granule. 

For  the  case  of  an  aluminum  or  magnesium  granule,  k  is  greater  than  k  (j.ss 
by  a  factor  of  1000.  Since  o'V'k  is  of  the  order  of  microseconds,  the 
time  required  for  the  establishment  of  thermal  equilibrium  is  about  5  billi¬ 
seconds.  By  a  similar  analysis,  we  find  the  time  required  for  the  estab¬ 
lishment  of  mechanical  equilibrium  is  of  the  same  order  of  magnitude. 

The  last  calculation  shows  that  if  the  smallest  time  scale  of  interest 
is  greater  than  10  milliseconds,  the  mixture  can  in  fact  be  treated  as  a 
system  in  thermodynamic  equilibrium.  On  the  other  hand,  if  the  smallest 
time  scale  is  of  the  order  of  1  millisecond,  we  may  define  a  mean  temper¬ 
ature  for  the  granule  and  another  for  the  gas.  The  same  applies  to  the 
velocity  of  the  granule  and  the  gas. 

To  be  sure,  we  may  visualize  a  situation  in  which  the  temperature 
and/or  velocity  in  the  gas  surrounding  the  granules  may  vary  considerably 
throughout  most  of  the  region  enclosed  by  the  small  volume  S?  and, 
in  spite  of  this,  a  mean  temperature  and  velocity  may  be  defined  by  some 
averaging  process  vithin  the  volume  .  If  such  is  the  case,  the 

smallest  macroscopic  time  of  interest  may  be  considerably  smaller  than  a 
millisecond.  However,  the  exact  nature  of  the  approximation  used  becomes 
rather  unclear. 

3.  Conservation  Laws  The  conservation  laws  governing  the  motion  and 
deformation  of  the  medium  can  be  easily  derived.  For  the  sake  of  simplicity, 
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we  shall  consider  the  granules  as  all  of  one  species  and  the  gas  as  another. 
Generalization  of  the  analysis  to  mixtures  of  many  species  does  not  present 
any  special  difficulties.  We  shall  use  the  superscripts,  the  "primes"  and 
"double -primes",  to  signify  quantities  pertaining  to  the  granules  and  the 
surrounding  gas  respectively.  Thus,  ^ 1  and  denote  respectively  the 

'i  .A  j 

mass  of  the  granules  and  of  the  gas  per  unit  volume  of  the  mixture j  u 
and  u"  denote  respectively  the  velocities  of  the  granules  and  of  the  gas. 
The  law  of  mass  conservation  can  he  expressed,  simply  as: 


+  V.(^''u/+ ^"u")  =  o  (3*1) 

^  t 

— X  —A 

If  and  if"  denote  respectively  the  "body  force  acting  on  the 

granules  and  on  the  ga.s  per  unit  volume  of  the  mixture,  the  equation  of 
momentum  balance  is 


.  f  u“)  +  V.  (?u‘  u'  .f  u"u")tVp  =  f  +  f "  (3.2) 

^  t 

where  p  denotes  the  gas  pressure. 

In  deriving  the  last  equation,  we  have  neglected  a  term  which  may  be 
called  the  "pressure  contributed  by  the  granules".  We  recall  that  the 
quantity  which  is,  on  the  macroscopic  scale,  called  "pressure"  is,  in 
actuality,  the  average  value  of  a  component  of  the  momentum  flux  on  the 
microscopic  scale.  To  be  more  precise,  if  we  denote  the  velocity  of  a 
gas  molecule  relative  to  an  observer  moving  with  the  mean  velocity  u"  by 
"c"  and  if  m"  is  the  mass  of  the  molecule,  the  pressure  tensor  pljj  is 
defined  by  the  phase  average: 


p« =<“'**> 


where  /  signifies  an  average  over  the  velocity  space.  The  scalar 

\  '  u 

pressure  p  is  by  definition  -5-  .  In  exactly  the  same  manner,  a 

pressure  contributed  by  the  random  motion  of  the  granules  may  be  defined 


in  terms  of 


where  c'  is  the  velocity  of  a  granule  relative  to  an  observer  moving  with 
the  mean  velocity  "u'  ;  m(  denotes  the  mass  of  the  granules;  and/} 

<W/T' 

signifies  an  average  taken  over  the  phase  space  spanned  by  C',  and  the 
temperature  of  the  granules  X1  .  However,  unlike  the  case  of  gas  molecules 
C'  is  usually  small  compared  with  uy  so  that  the  "pressure"  p  s  -4-  is 

in  general  much  smaller  than  p*  u’2  •  For  this  reason,  we  shall  neglect  p 

the  "pressure"  due  to  the  random  motion  of  the  granules. 

Under  the  same  approximation,  we  may  write  down  the  energy  equation: 

=  (3.3) 


4.  Equations  Governing  the  Motion  of  the  Granules  and  the  Gas  Let  us 

/  If 

define  <0  and  «  by  the  formulas: 

\P'  _  ,-v  * 

JF  +  ^  u  =  w  (4.i) 

if+V.  p"u"=6)"  (4.2) 

H 


Evidently  w  has  the  significance  of  being  the  rate  of  increase  of  the 
total  mass  of  the  granules  per  unit  volume  of  the  mixture,  while  ti*U  is  the 
rate  of  increase  of  the  mass  Of  the  gas  per  unit  volume  of  the  mixture. 

Eqs.  (4,1)  and  (4.2)  thus  represent  the  equations  of  mass  balance  of  the 
two  components.  Furthermore,  by  virtue  of  (3«l)  that 

Oi-f  tu’ s  0  ('ll  o\ 


In  a  similar  manner,  we  define  and  F,(  by 

l£iVv.(P"  i «") 

b  t 


(4.4) 

(4.5) 
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V'  is  thus  the  total  force  everted  by  the  gas  on  the  granules  per  unit 
volume  of  the  mixture  and  F  is  the  reaction  of  this  force  on  the  gas. 
Evidently 

FVF"=o  (4.6) 


Furthermore,  eqs»  (4.4)  and  (4.5)  may  be  interpreted  as  the  equations  of  j. 
motion  of  the  granules  and  of  the  gas  respectively. 

Finally,  we  define  and  by  the  formulas: 

tfitfaH-'l  +v.[f'^a't  i/j.  (fVF').Sv?-  (t.7) 

If  we  interpret  (4.7)  and  (4.8)  as  the  energy  equations  for  the  granules 
and  the  gas,  the  quantity  c ^  will  then  have  the  significance  of  being 
the  heat  flow  into  the  granules  per  unit  volume  of  the  mixture,  while 
is  that  flow  into  the  gas  per  unit  volume  of  the  mixture.  From  (3.3), 
we  see  that 


r-'J  -fill 

t-  r  •  u 


O 


(4.9) 


Eqs.  (4.3),  (4.6),  (4.9)  provide  three  relations  for  the  six  variables: 

F'  ,  f"  *  We  need  three  more  relations  before 

the  problem  is  completely  determined.  In  general,  these  relations  must 
be  postulated.  However,  when  the  state  of  the  mixture  departs  only  slightly 
from  the  state  of  local  thermodynamic  equilibrium,  the  form  of  those  re¬ 
lations  can  be  determined  from  thermodynamic  considerations.  (See  the  next 
section.) 

We  note  that  eq.  (4.4)  can  also  be  written  as: 


r  at* 


+  Ck)  ^ 


where 


d_ 

at1 


L. 
1 1 


rj 

u  .  V 


Likewise,  if  j  .  —  ^ _  +  u,,»  V 

at"  ~ bt 


(4.10) 


we 
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have: 


oh*£!  +  ^"u'-vp  r  +-  fh 

v  j  t" 

Similarly  we  can  rewrite  (4.7)  and  (4.8)  as 


(4.n) 


?,L(nViu,1)1.w|(n+i:>A'2j=(|VF').M,+  1'  .  (4.i2) 

d  t 

( n'V 1 U»> -  V.(p5")+  (|?VF»).u"+<j"  (4.13) 

Furthermore,  if  we  eliminate  fVF'  from  (4.10)  and  (4.12),  we  obtain 

(4.14) 

d  t 

Likewise,  from  (4.1l)  and  (4.13) 

f«4S;tw"(TT-r  u"1)=-pv.t»+«{"  (4.i5) 

a  l 

The  last  equation  may  also  he  rewritten  as 

?“[$'  *  f  )+»*(«"-  f  u“2)=  ?"  ^ 


(4.12) 


(4.15) 


(4.16) 


where  use  has  been  made  of  (4.2).  In  terms  of  the  entropies,  eqs.  (4.l4) 
and  (4.16)  assume  the  simple  form 


(4.17) 

i 

(4.18) 

since 


T'J5'=<1U’  ,  a’=H' 

for  a  rigid  granule  and  T-JStJUtpJ(^  for  a  gas. 


(4.19) 


H  m” 


5.  Entropy  Production  Let  ^  be  the  density  of  the  miscture.  Then 

F-  S.  A  Ell?  OEDBABOE 
AiJEEDEEi?  PECVIIG  GEOUJJD,  Mil 
-9-  o&DBfi-n, 


It  is  convenient  to  introduce  the  mass  fraction  y '  which  is  defined  as 
the  total  mass  of  the  granules  in  a  unit  mass  of  the  mixture.  It  follows 
that 

(5.2) 

Likewise  if  y"  denotes  the  total  mass  of  gas  in  a  unit  mass  of  the  mixture, 
we  have 

Y"=e"/?  (5-3> 


Obviously, 


Y,+  Nf"  =  l  (5.4) 

We  shall  also  define  the  velocity  u  of  the  mixture  by  the  formula 

■  -N  I  k  |  \l  — ^  1 1  f  i 

^  u  o\'  (5.5) 


Eq.  (4.1)  can  then  be  rewritten  as 


|  y .  t  ^  =  o 


(5.6) 


which  is  the  continuity  equation  for  the  mixture.  Finally,  let  us  in¬ 
troduce 


ot  n 


(5.7) 


and 


v'=  u' - u 


-A  II 

V 


-  u  _  u 


(5.8) 


It  may  be  easily  verified  that  (4.1)  and  (4.2)  may  be  rewritten  as: 

214av.(pY,''')  =  («,/p 

U  u  \ 

eiV-jrV.tfY'v")=“7e 

Dt  \  1 
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Also  (4.18),  (4.19)  136001116: 

r$*rAvs=^[VU7«HJ“>] 

If  5  denotes  the  entropy  of  a.  unit  mass  of  the  mixture, 

5  =  Y'5'  +  Y"S" 


(5*11) 

(5*12) 


(5.13) 


and  we  may  compute  £5  from  (5.9),  (5.10),  (5.1l),  (5.12)  and  (5. 13) 

ot 

Omitting  the  algebraic  details,  we  find  ultimately  the  equation 


?a|+v.(^'S'tfllv«S")=  § 


(5.14) 


where 


(5.15) 


y/  1  y/f  1  yi  v  2  '  y 

andG-*  ,  G-'1  are  Gibbs’  free  energy  per  unit  mass  of  the  granules  and 
of  the  gas  respectively.  Eq.  (5.14)  has  a  simple  meaning,  for  it  may  he 
rewritten  as 


eu5"S“)=  § 


(5*16) 


where  use  has  been  made  of  (5*6)  and  (5.8).  If  we  integrate  (5.16)  over 
a  region  enclosed  by  a  smooth  control  surface,  we  have 

d  <3 


pf  Jv  *  fet'S' + y  s  *).a/\ = ^ 


(5.17) 


The  first  term  of  this  equation  is  clearly  the  rate  of  increase  of  entropy 
inside  the  volume.  The  second  term  represents  the  amount  of  entropy 
convected  out  of  the  region.  The  last  term  must  then  be  the  entropy 
produced  in  the  region.  Consequently,  signifies  the  entropy  produced 
in  a  unit  volume  of  the  mixture.  Making  use  of  (4.3),  (4.6)  and  (4.9), 


we  have: 


a  - 

6.  Small  Departure  from  Thermodynamic  Equilibrium  At  thermodynamic 


equilibrium. 


u*  =  u" 


F'  -=  F*«  O 


T'=  T" 


<r=&" 


V  =  f  =  ° 


Gl'  —  40U—  O 


(6.1) 


so  that  ^=0  and  the  entropy  of  the  system  attains  a  maximum.  When  the 
system  is  near  equilibrium,  we  can  expect  the  following  linear  relation 
between  the  "fluxes"  and  the  "forces". 


F'=-F"  =  k 


(6.2) 


i  =  -f-  (6.31 

\  /  1. 

where  X)0,  |4>0,  h^O  andA=  A-  particular,  if  there  is  no  reaction 

between  the  granules  and  the  surrounding  gas,  we  have  simply: 


(6.4) 


F,=  -F"  =  k(u"-uO 
1,si"  =  V(T'-T') 


(6.5) 


On  the  other  hand,  if  the  smallest  macroscopic  scale  is  large  compared 
with  the  time  required  for  the  granules  and  their  surrounding  gas  to  attain 
mechanical  and  thermal  equilibrium,  then  u"  and  Tf  =  Tn  .  For  this 

case,  if  there  is  mass  exchange  between  the  granules  and  the  gas  surrounding 
them,  we  have 


=  -  w"  =  -if.  (&"-&') 


(6.6) 


Of  course,  in  nearly  all  cases  where  chemical  reaction  enters  as  an  im¬ 
portant  factor,  G-11-  &■  is  not  small  and  the  system  undergoes  changes 
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which  are  very  far  from  thermodynamic  equilibrium.  In  such  cases  (6.6) 
cannot  be  used  and  must  be  replaced  by  an  appropriate  rate  law  determined 
from  other  considerations c  On  the  other  hand,  the  relaxation  times  between 
the  granules  and  the  surrounding  gas  are  in  general  small  and  the  state 
of  the  system  is  never  too  far  from  the  state  of  thermodynamic  equilibrium. 

This  means  that  (6.2)  and  (6.3)  may  be  used  in  most  practical  applications. 

7.  Characteristic  Theory  The  basic  equations  governing  the  one-dimensional 
motion  of  the;  mixture  are  summarized  below: 


a!  -  *>• 

at  a* 

(7»l) 

at  <)* 

(7.2) 

H'Oevd  f<+F/ 

(7-3) 

SfV  y>p  f.+  p. 

at  a*  a*  * 

(7.4) 

1  u'JS’  _  $y?' 
st  it  - 

(7.5) 

St  S*  ~  ^  '  ' 

(7.6) 

where 

(7.7) 

(7.8) 

In  addition,  we  have  the  equation  of  state. 

This  can  be  readily  shown  to 

"  RT" 

(7.9) 

where  R  is  the  gas  constant  and  T  is  the  temperature  of  the  gas.  Further¬ 
more,  the  entropies  of  the  granules  and  the  gas  are  respectively 
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(see  ref.  2)  To  calculate  the  characteristics  of  the  system,  we  shall 


denote  the  change  of  any  variable  along  a  characteristics  by 

+CL 

St  at 


«§  t 


Thus 


(7.12) 


where  c  denotes  the  speed  of  propagation  of  the  characteristics  surface. 
Let  us  now  eliminate  the  spatial  derivatives  from  all  of  the  equations 
and  look  for  the  specific  directions  along  which  the  variation  of  the 
different  quantities  can  be  determined  by  an  equation  in  this  direction. 
These  are  the  characteristics  directions.  Thus,  from  (7«5)j  we  find  the 
characteristics  direction  in  the  xt  -space  is  specified  by  the  wave  speed 
C=  U'  and  along  which  =$'/?'  •  Carrying  out  this  procedure,  we 
obtain  the  following  results.  There  are  six  sets  of  characteristics: - 

l)  A  triplet  of  the;  characteristics  Csia*  along  which 
) 


jg«  ^ 


(7.13) 


and 


=  $VF'-  u'lo1  (double)  (7.1*0 

2)  A  single  family  of  characteristics  propagating  at  c=  u"  along  which 

If  =  r 1  ?*  (7.i5) 

3)  Two  families  of  characteristics  propagating  with  the  speeds 
along  which  we  have 


St 


_ fc  $7frt-uBSS!l  = 


o 


(7.16) 


Because  of  the  degeneracy  of  the  characteristics  ralations  (7.1*0 j  we  have 
only  five  equations  for  the  six  variables : 5  ,  b" ,  \a'  ,  ,  p*  and  p" 

Note  that  the  pressure  can  always  be  calculated  from  the  equation  of  state 
and  expressed  in  terms  of  these  variables. 


To  complete  the  system  we  need  an  additional  equation.  This  is 
provided  by  (7.1)  if  we  eliminate  jL_  from  it  with  the  help  of  (7. 12). 


Thus  we  find  that  along  c  =  u 


£I'+d'Iu,=  w'  (7.17) 

St  x  )x 

These  relations  are  adequate  for  studying  wave  problems  in  the  mixture  by 
numerical  methods. 


8.  Small  Amplitude  Waves  For  small  amplitude  waves  we  may  linearize 
the  governing  equations  and  an  analytical  solution  may  be  constructed.  In 
this  calculation,  we  shall  assume  oi'=  0  .  The  linearized  form  of  the 
governing  equations  are; 


+  p'  V.  u's  O 

it  V 

—  +  ?!  u"=  o 

}  t  u 

+0"h''-  _7P 

r°<)t 

p'r'  iV  .  p "c"£L"-^p 

-  ?.cPJt  -it 

X(T--r) 

£  p"  tv 

p.  ~m r  t: 

where  subscript  ^o11  characterizes  the  undisturbed  field.  We  have 
assumed  that  the  body' forces:  =  ss  .y  0  . 

this  section  p  y  u‘  t  u",  p1 ;  p1'  ,  T '  f  T"  et:c?^r^,'| 


(8.1) 


(8.2) 


(8.3) 


(8.4) 


(8.5) 


(8.6) 


(8.7) 


.In 

iJBIYSHEBB  PBOTOIG  ©BOMB.  &»• 
OBPBGhTt 


to  be  the  perturbed  quantities, 


If  X»' ,  J'  ,  p“  and  T "  are  eliminated  from  the  above  system  of 


equations,  we  have 


ht 


where 


°<  =  k(^  +  yr) 

=  k  I 

^ x  +  Tu7^ 


4  =  >( 


hr) 


av  ?icv 


Eliminating  u  we  obtain 


where 


at1 


Cl=  I /of 

t,=  i/5 


a.=  aJ&£ 


a2=a„ 


«+  "S 

'  *5 


(8.8) 


(8.9) 


(8.10) 


(8.11) 


(8.12) 


(8.13) 


(8.14) 


(8.15) 


(8.16) 


(8.17) 


(8.18) 


(8.19) 


and  ao,\e f”  '  sound-  speed  in  the  gas  when  there  is  no  granules. 

The  solution  of  (8.14),,  satisfying  the  initial  condition  that  the  medium 


is  undisturbed  at  t-  o  and  the  boundary  condition  associated  with  the 
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piston  problem,  i.e,; 


u‘=  u.  H(t) 


(8.20) 


at  XiO 


r  0  TI  L 


I  «+  s 


£  t-A* 


2TU  J  A 


(8.21) 


where 


A  =  — 

a°f(s+P>)(s+  T,) 


(8.22) 


For  large  S  , 


s  t-Ax  =  s(  t-  —)  - 


(°<  -  ft)  +  (5-%)  x  r>  (  '  ) 

2  30  '  \f» 


(8.23) 


Hence,  if  x)  aQt  the  integral  must  be  closed  to  the  right  of  the  imaginary- 

axis  and  hence,  o  .  For  x  <  aot  ,  the  integral  must  be  closed  to  the 

left  and  f>  is,  in  general,  different  from  zero.  In  fact,  the  value 
of  as  x  — >  a Ot-o  is  easily  computed  by  applying  a  Tauberian  theorem. 


Thus,  it  is  found  that 


Uvn  JP_  expLisld^LtSiSi-tl 

X-9a„t-o*p0  Ba  ’  L  *  ■* 

—  Ms  f-  1-  ( il-  A—  \  tj 

a<>  »  p;c'c 1  i 


(8.24) 


Hence  the  wave  front  is  rapidly  damped.  This  is  as  it  should  be;  for, 
the  wave  front  represents  a  signal  propagating  between  the  granules.  The 
probability  that  it  is  not  scattered  by  the  granules  decreases  exponentially 
with  time. 

To  examine  the  behavior  of  p  at  the  piston  for  large  t  ,  we 
introduce  a  new  variable  S-=ts  .  Thus,  at  the  piston. 


_  f0  u0ao  1  I s'+  ^ s'+t,  t '  gb 

2TTI  K  s'+-"$  t 

It  CO 


(8.25) 


so  that  as  t— *  oo 


fc> 

1 


pM  u  a0  /_2L  -SJ 
V(b  1 


(8.26) 


Hence 


Urn  (  \  z  Up~  (8“27) 

if  MU  Pe  &i 

where 


v  _  y  £i  (8.28) 

0  Hu  "  3  ^ 

The  last  formula  shows  that  for  large  time  the  pressure  level  is  determined 
by  the  equilibrium  sound  speed  a,  .On  the  other  hand;  for  small  time 
the  pressure  at  the  piston  is  determined  by 


tovi  _  _U*  f  \  _  (ot-ft)  t  4-  O C  t2-)!  (8.29) 

X  ft,  ~  Bo  l  1  -* 

Finally;  the  wave  profile  at  large  time  can  also  be  deduced  by  using 

the  method  of  the  steepest  descent.  (For  details;  see  ref.  5) • 

Let  us  introduce  the  function 


W(5l=  s[ 


then 


w 

'  2TT  i 


LOT) 


A. 

'.Loo 


_x_  1 

[  s+*  s-*-V 

*ot  V 

S+  S+  tj 

t  V/Cs) 

o(vS 

e  d 

(8.30) 


(8.31) 


As  t  — *  00  ;  the  main  part  of  the  contribution  of  the  integrand  is  con¬ 
centrated  near  the  point  where  W(s)  is  a.  maximum.  Since  v/  is  an 
analytic  function  of  S>  }  this  point  is  in  actuality  a  saddle  point.  A 
careful  examination  of  the  trajectory  of  the  relevant  saddle  point  shows 
that  it  is  located  on  the  real  axis  and  crosses  S  =  D  at 


X  _ 

a„t  “  V 


a* 

x  "  ae 


(8.32) 


Now  the  point  S-  o  is  a  pole  of  the  integrand;  it  follows  that;  for  *<a2t 
the  solution  is  of  one  type  while;  for  X>a2t  f  the  solution  is  of  a 
different  type.  Taking  into  account  the  confluence  of  the  saddle  point 
and  the  simple  pole  at  5=0  ;  we  find 
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(8.33) 


t*  -  j_  Up  \S01 

'3,|  5. 

where  is  the  saddle  point  and  Since  Se^  o  for  X<3ot 

and  SeK.o  for  x  }  aat  ,  the  wave  profile  consists  of  a  small  precursor 
wave  propagating  at  the  sound  speed  a„  which  is  followed  by  the  main 
signal  propagating  at  the  equilibrium  speed  .  This  conclusion  is 

qualitatively  similar  to  the  result  of  G.  F.  Carrier  (ref.  8)  who  considered 
the  same  problem  without  thermal  relaxation.  The  main  difference  is  that 
the  main  signal  propagates  with  the  equilibrium  sound  speed  adhere.  We 
note  that  the  same  conclusion  was  deduced  earlier  in  ref.  2  where  a 
Thermodynamic  analysis  was  made. 


-19- 


REFERENCES 


1.  Hershkowitz,  J.,  Schwartz  F.  and  J.  V.  R.  Kaufman.  "Combustion  in 
Loose  Granular  Mixtures  of  Potassium  Perchlorate",  Eighth  Symp.  (int.) 
on  Combustion,  i960. 

2.  Chu,  B.  T.,  "Thermodynamics  of  a  Dusty  Gas  and  Its  Application  to 
Some  Aspect  of  Wave  Propagation  in  the  Gas".  Brown  University  Rept. 
DA476l/l  United  States  Army  Ordnance  Corps,  Picatinny  Research  Center 
Cont.  No.  DA-19-020-0RD-4761.  (i960) 

3.  Chu,  B.  T.,  "Thermodynamics  of  a  Porous  Elastic  Solid  Containing  a 
Compressible  Fluid".  Brown  University  Rept.  DA476l/2,  United  States 
Army  Ordnance  Corp.,  Picatinny  Research  Center,  Cont.  No.  DA-19-020- 
0RD-4761.  (i960) 

4.  Chu,  B.  T.,  "Thermodynamics  of  a  Loosely-Packed  Granular  Aggregate", 

Brown  University  Rept .  DA476l/3j  United  States  Army  Ordnance  Corp., 
Picatinny  Research  Center,  Contract  No.  DA-19-020-0RD-4761.  (1962) 

5.  F.  A.  Williams,  "On  the  Statistical  Description  of  Sprays".  Phys.  of 
Fluids,  1958,  1,  541. 

6.  F.  A.  Williams,  "Spray  Combustion  Theory"  Eighth  Symposium  on  Combustion 
(international),  Williams  and  Wilkins,  1961. 

7.  F.  A.  Williams,  "Detonation  in  Dilute  Spray",  in  Detonations  and  Two- 
Phase  Flow.  A.R.S.  Progress  Series  in  Astronautics  and  Rocketry.  Academic 
Press,  Inc.,  N.  Y.,  1961. 

8.  Carrier,  G.  F.,  "Transmission  of  Blast  Through  a  Dust-Filled  Gas — 
Linearized  Treatment",  BRL  Memo  Rept.  #873*  Aberdeen  Proving  Ground, 
Maryland. 


-20- 


